Also noting here this relevant research article Synthetic lethality in large-scale integrated metabolic and regulatory network models of human cells.

Also noting here this relevant research article Synthetic lethality in large-scale integrated metabolic and regulatory network models of human cells.

yes, it mentioned another method gMCStool that can predict essential genes in context-specific models.

Other checks than essentiality are also welcome

@mihai-sysbio Could you explain the reason for the failing task?

@mihai-sysbio Could you explain the reason for the failing task?

The task isn't actually failing, really. As mentioned in the comment above, it's simply that the output of the gene essentiality is so long that it cannot be posted as a PR comment. I couldn't figure out how @haowang-bioinfo ran the calculations to obtain a much more compressed output. As visible in the workflow changes, the code here is just calling functions defined elsewhere.

I couldn't figure out how @haowang-bioinfo ran the calculations to obtain a much more compressed output. As visible in the workflow changes, the code here is just calling functions defined elsewhere.

the concise output was an outcome of a series of manual operations of running functions, copy/paste results, and modifications with text editor

If these last manual steps could be ran as a function, then we would be ready to run the gene essentiality analysis on every PR 💪🏻

As is suggested for MACAW results in #840, a summary output could be automatically posted in the PR, while a more detailed text file could be written and committed as part of each PR. This could be organized in the same folder as the MACAW results, and a README.md file could indicate for both output files from which PR they derived from.

This PR has been automatically tested with GH Actions. Here is the output of the MACAW test:

Starting dead-end test...
 - Found 1514 dead-end metabolites.
 - Found 1319 reactions incapable of sustaining steady-state fluxes in either direction due to these dead-ends.
 - Found 1977 reversible reactions that can only carry steady-state fluxes in a single direction due to dead-ends.
Starting duplicate test...
 - Skipping redox duplicates because no redox_pairs and/or proton_ids were provided.
 - Found 447 reactions that were some type of duplicate:
   - 0 were completely identical to at least one other reaction.
   - 13 involve the same metabolites but go in the opposite direction or have the opposite reversibility as at least one other reaction.
   - 447 involve the same metabolites but with different coefficients as at least one other reaction.

This and a more detailed output from MACAW are also committed to data/macawResults/.

Note: In the case of multiple test runs, this post will be edited.

This PR has been automatically tested with GH Actions. Here is the output of the gene essentiality test:

     cellLine     TP     TN     FP    FN     accuracy    sensitivity    specificity       F1         MCC   
    __________    __    ____    __    ___    ________    ___________    ___________    ________    ________

    {'DLD1'  }    36    2185    59    279    0.86792       0.11429        0.97371       0.17561     0.15291
    {'GBM'   }    34    2165    61    298    0.85966       0.10241         0.9726       0.15925      0.1333
    {'HCT116'}    46    2207    53    309    0.86157       0.12958        0.97655       0.20264     0.19047
    {'HELA'  }    30    2263    69    254    0.87653       0.10563        0.97041       0.15666     0.12398
    {'RPE1'  }    14    2204    81    259    0.86708      0.051282        0.96455      0.076087    0.025853
    {'all'   }     7    2408    92    109    0.92317      0.060345         0.9632      0.065116      0.0254

Note: In the case of multiple test runs, this post will be edited.

Finally 😄 it seems to work! Ready for review.

This is amazing 🤩 thank you @edkerk for seeing this through 💪🏻